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Abstract 

We propose a new one-parameter family of Landau gauges for Yang-Mills theories which can be formulated by 
means of functional integral methods and are thus well suited for analytic calculations, but which are free of Gribov 
ambiguities and avoid the Neuberger zero problem of the standard Fadeev-Popov construction. The resulting gauge- 
fixed theory is perturbatively renormalizable in four dimensions and, for what concerns the calculation of ghost and 
gauge field correlators, it reduces to a massive extension of the Fadeev-Popov action. We study the renormalization 
group flow of this theory at one-loop and show that it has no Landau pole in the infrared for some - including 
physically relevant - range of values of the renormalized parameters. 
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1. Introduction 

In his seminal work on non-abelian gauge theories 
0], Gribov showed that the issue of fixing a gauge is far 
more complicated than in the abelian case and may be 
hampered by some ambiguities. For instance, for SU(A0 
Yang-Mills (YM) theories, the Landau gauge conditionj 
dfiAfj = is satisfied by many field configurations that 
are equivalent up to gauge transformations. To com- 
pletely fix the gauge one has to specify how to deal with 
these Gribov copies. 

For Landau gauges, Gribov copies are given by the 
extrema of the functional F[A,U] = J d d xtr(A u ) 2 with 
respect to U for a given field configuration A, where 
A u is the transform of A under the gauge transforma- 
tion U. The standard Faddeev-Popov (FP) procedure 
corresponds to summing over all extrema of F[A, U] 
(minima, maxima and saddle points) with a minus sign 
for the Gribov copies whose FP operator has an odd 
number of unstable directions QJJ] . For compact gauge 
groups, such as SU(A0, this alternation of sign leads to 



'We consider the Euclidean theory in d dimensions. When neces- 
sary, a lattice discretization is implicitly considered. 



an exact compensation of contributions from the vari- 
ous Gribov copies for gauge invariant quantities. As 
a consequence, formal expressions of physical observ- 
ables in fact appear as an undetermined 0/0 quotient f^j]. 
This so-called Neuberger zero problem forbids the use 
of e.g. lattice regularizations to devise a nonperturba- 
tive version of the gauge-fixed theory. It is usually dis- 
regarded in the high-energy, perturbative regime, where 
Gribov copies are thought to be harmless and are sim- 
ply ignored. Another drawback of the FP gauge-fixing 
is that the resulting perturbation theory shows a Landau 
pole in the infrared (IR) - where the running coupling 
constant diverges - and is thus useless for studying the 
low energy properties of the theory. 

Several other Landau gauge-fixing schemes can be 
considered. One can for instance pick up a particu- 
lar minimum of F [A, U] for each field configuration A. 
This has the advantage of being relatively easy to im- 
plement for numerical calculations on the lattice J^. 
However, depending on the algorithm used to find the 
minimum, different Gribov copies will be attained, cor- 
responding to different choices of gauges. This renders 
the comparison between different calculations (of non 
gauge invariant quantities) somewhat tricky. In partic- 
ular, it has been shown that different choices of algo- 
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rithms may lead to a variation of the ghost propagator 
by up to 10% at low momentum fl. 

One way to cope with this issue is to choose the ab- 
solute minimum of F[A, U] for each A, the so-called 
absolute Landau gauge 051]. This is clearly a pretty hard 
numerical task since this functional typically presents 
numerous local minima. Gribov and Zwanziger [l], 0] 
proposed to restrict the functional integration over the 
A field to the first Gribov region (corresponding to the 
minima of F[A, U]). However, this region is not free 
of Gribov ambiguity [7]. Recently, a refined version of 
this Gribov-Zwanziger approach was shown to lead to 
predictions for the gauge field correlation functions that 
agree well with lattice simulations [8]. 

In this letter, we propose a new one-parameter fam- 
ily of Landau gauges which are free of Gribov ambi- 
guities and which turn out to be simple to implement 
for practical perturbative calculations. These are good 
gauge-fixings in the sense that gauge-invariant quanti- 
ties are independent of the gauge fixing procedure. We 
also point out that this family of gauges interpolates be- 
tween the FP and absolute Landau gauges. Another one- 
parameter family of Landau gauges was proposed in 01 . 

One important aspect of the gauge fixing procedure 
proposed here is that it can be fully implemented by 
means of standard functional integral techniques and is 
thus well suited for analytic calculations. Indeed, our 
procedure amounts to taking a particular average over 
Gribov copies of each gauge field configuration. This 
average can be dealt with by means of a standard replica 
trick - borrowed from the theory of disordered systems 
in condensed matter physics - which, together with a 
Fadeev-Popov (FP) construction, allows for a simple 
path integral formulation. Another key aspect is that our 
procedure does not suffer from the Neuberger zero prob- 
lem of the standard FP construction. It is in this sense 
a fully justified field-theoretical description of YM the- 
ory which takes into account the Gribov ambiguity from 
first principles. 

We show that the model is perturbatively renormaliz- 
able in four dimensions thanks to the underlying sym- 
metries. We also demonstrate that, within perturbation 
theory, correlation functions of YM fields are identical 
to those obtained in a massive extension of the usual 
FP gauge-fixed action, which is a particular case of the 
Curcci-Ferrari (CF) model ifToh . This provides a solid 
first-principle field theoretical justification for the phe- 
nomenological approach recently proposed in 1 1 ill 1 211 . 



One-loop calculations in this model were shown to re- 
produce quantitatively lattice results for the gluon and 
ghost propagators in the Landau gauge. Finally, we 
discuss the renormalization group flow of the present 



gauge-fixed YM theory at one-loop and show that it is 
free of Landau pole for some range of renormalized pa- 
rameters. 

2. The gauge-fixing procedure 

The main idea of this work is to deal with Gribov 
copies not by choosing a unique one as is often done, but 
by averaging over them with a given weight. To com- 
pute usual YM correlators, we first average over Gribov 
copies for each given gauge field configuration A and 
then perform an average over the gauge field configura- 
tions with the usual YM weight. 

To be specific, we consider a SU(N) YM theory in 
the Landau gauge. The classical action reads 



Sym[A] 



= ^J}r(;v)\ (i) 



where = d^Ay-dyA^-igolA^Ay], JT = J d d xandg 
is the (bare) coupling constant. Here and below, when 
a field is written without color index, a contraction with 
the generators of the group in the appropriate represen- 
tation is understood. Our convention for the generators 
is tr T a T b = \5 ab and [T a , T h ] = if abc T c . 

For any operator 0[A], we define the average over the 
Gribov copies of a given field configuration A as: 



= ^OlA^e-^ 

2, sO>- Sw[A ' ] 



(2) 



with weight factor 

Sw[A] =A>J\r(A / ,) 2 . (3) 
and where the gauge transform A u is defined as: 



UA^U 



~* + —UdM- 1 . 

go 



(4) 



In Eq. ©, the sum runs over all Gribov copies, that is 
over all t/, which satisfy d^A^ = or, equivalently, 
over all extrema of Sy/[A Ui ], for a given A; s(i) is the 
sign of the functional determinant of the FP operator 

{5 ab d 2 -gof abc A^)6 (d \x-y) (5) 

taken at A - A Ui . Finally, fio is a free parameter. The av- 
eraging procedure (fJJ-OJ is inspired from the one pro- 
posed in II 1 3fl in the context of the Random Field Ising 
Model. Note that a similar average over Gribov copies 
was considered in lfl4ll to implement a simulated anneal- 
ing as an efficient method to approach the absolute Lan- 
dau gauge. Here, however, the sum runs over all Gribov 
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copies, not only on those in the first Gribov region. No- 
tice though that, for fio not too small, the copies outside 
the first region are suppressed by the weight factor (0. 

Actually, the limit — > oo selects the absolute mini- 
mum U - f/abs of Sy/[A U ] and our averaging procedure 
simply corresponds to the absolute Landau gauge: 



lim (0[A]> = 0[A V "*] . 



(6) 



In the opposite limit fio — > 0, all Gribov copies con- 
tribute the same in the average (0 up to the sign fac- 
tor s(i). Since there are as many contribution with each 
sign, the denominator in vanishes: ■ s(i) = 0. For 
any fio > 0, the degeneracy is lifted and the denomina- 
tor in (|2j differs from zero in general^ which solves the 
Neuberger zero problem. 

Once the average over Gribov copies has been 
performed for each individual field configuration we 
perform the usual average over YM field configurations, 
hereafter denoted by an overall bar: 



0[A] = 



J DAO[A]e~ s ™ [A] 



,-S ymL^] 



(7) 



To summarize, computing a given observable O in our 
gauge fixing procedure amounts to first average over 
Gribov copies and then to average over YM field con- 
figurations, that is 

<OA]>. (8) 

A crucial remark is in order here: observe that gauge- 
invariant operators such that 0- mv [A u ] = (9; nv [A], are 
blind to the average (0: (<9j nv [A]) = <9i nv [A], which 
guarantees that our gauge fixing procedure does not af- 
fect physical observables. In particular, one has 



(O im [A])=O im [A]. 



(9) 



It is crucial to introduce the denominator in 01 in order 
for this fundamental property to hold. 

3. Functional integral formulation 

The notations used in the previous section where cho- 
sen to make contact with the physics of disordered sys- 
tems 1 15], the theoretical description of which also typ- 
ically involves a two step averaging procedure: One has 
to average first over thermal fluctuations for a given re- 
alization of the (quenched) disorder and then over the 



"Strictly speaking, we cannot exclude that for a given value of /Jo 
some field configurations yield a vanishing denominator. We assume 
such field configurations to be of zero measure. 



disorder. In the above proposal for YM theories, the 
analog of the first (thermal) average is the one over Gri- 
bov copies, Eq. (0, for fixed gauge field A. The latter 
thus plays the role of the disorder field to be averaged 
over in the second step. Although such two-step aver- 
ages are rather easily implemented in numerical simula- 
tions, they are technically challenging for analytical cal- 
culations. The method of replica is designed to handle 
this issue. We show here how to adapt it to the problem 
at hand. 

Let us first rewrite the average over Gribov copies as 
a functional integral. Following the standard FP proce- 
dure, the sums over Gribov copies appearing in 01 - 
including the sign s(i) - can be represented by a func- 
tional integral over a SU(A0 matrix field U, FP ghost 
(Grassmann) fields c and c and a Lagrange multiplier h. 
One has, in particular, 



■S w [A u ,c,c,h] 



(10) 



where V(U, c, c, h) = DUDcDcDh, with DU the Haar 
measure on the gauge group, and 

S FP [A, c, c,h]=f {8^(3,0° + gof hc Ay) + ih a d u A;} 

(11) 

is the usual FP action. Collecting the set of FP fields U, 
c, c and h in a single symbol *V - we shall see shortly 
how this can be realized explicitly in a supersymmetric 
formulation -, the average (0 over Gribov copies can 
be written as 



(0[A]) 



J lXVO[A u ]e- s ™ [A 



j irvt 

with the gauge fixing action 

S C f[A,V\ =Syf[A u ] +S w [A u ,c,c,h]. 



(12) 



(13) 



The denominator of the previous expression can be 
formally rewritten using the replica trick 11511 



1 



.-Sgf[A,T] 



J k=\ 



(14) 

Here, the limit is to be understood as the value of the 
(analytically continued) function of n on the right-hand- 
side when n — > 0. The average over the disorder field A 
can then be formally written as 



(0[A]) = lim 
«->o 



■s[A,m] 



f*>Ai 



(15) 
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where 



S[A,{<V}]=S Y M[A]+Y J S GI! [A,<V k ]. (16) 



i<=\ 



Finally, using Eq. ( fT5l ) with 0[A] = 1, one can rewrite 
( fT5b in the more convenient form 



where the factor g l/2 (8, 6) = /3 O 60 - 1 on the right-hand- 
side can be seen as the invariant measure associated with 
curved Grassman dimensions, as discussed in Ref. {l(h . 

To gain more insight about the structure of the present 
theory, we write, alternatively, 



{0[A]) = lim 



—° j»A(nLi^)«" S[ " 



mi 



(17) 

Here, the choice of the replica k — 1 is arbitrary because 
of the obvious symmetry between replicas. 

It is useful, for perturbative calculations, to explicitly 
factor out the volume of the gauge group j DU. This 
can be done e.g. by performing the change of variables 
A -> A Ul and U k -> U k U~ l , V k > 1 in <[T7). Renaming 
(c\,c\,h\) — > (c,c,h), we get 

•S[A,c,cA{T)] 



<0[A]> = lim 



,c,c,h,{'V]] 



J D(A,c, c,h,{<V})0[A] e 

«" o J ©(A, c , c, ft, {-y}) [A 

(18) 

with £>(A, c, c, ft, {<V}) = £>(A, c, c, ft) (n* =2 and 
5 [A, c, c, ft, f/V}] = S ym[A] + 5 W [A] + S FP [A, c, c, ft] 

n 

+ Y J {Sw[A u "]+S w [A u \c k ,c k ,h k \). (19) 

As mentioned previously, the present theory can be 
given an elegant supersymmetric formulation, which 
proves powerful. To this purpose, we specify the sym- 
bol *V introduced before as the following matrix super- 
field: 



(23) 

where we introduced the covariant derivative D/V = 
d/V + igo'VAfi. We see that the action (O, with (Qj) 
and (|2TI) -(|231). thus describes a collection of n gauged 
supersymmetric nonlinear sigma models jvh . It is in- 
variant under the super gauge transformation A — > A w 
and 'V/t — > r V k r U 1 , ~ik — 1, . . . , n. This symmetry gets 
explicitly broken after one replica is singled out to ex- 
tract the volume of the gauge group. The action ( fT9l 
describes a gauged-fixed theory for n — 1 gauged su- 
persymmetric nonlinear sigma models with CF gauge 
fixing 5w[A] + Spp[A, c, c, ft]: 

S[A,c,c,h,{ ( V\\ = S Y m[A] +S w [A] + S FP [A,c,c,h] 



j_, Jx,r 



So 



(24) 



k=2 



4. Renormalizability 

We now prove the perturbative renormalizability of 
the model (fT9l)- (l24l) in d — 4. This is a rather non-trivial 
result given the presence of nonlinear sigma models, 
which are in general renormalizable in d — 2. For sim- 
plicity, we focus on the SU(2) case for which T a = cr°/2 



<V(x, 6, 6) = exp [ig (Sc + c9 + M)} U (20) with cr" the Pauli matrices and f ahc = e abc the Levy- 



living on a superspace made of the original d- 
dimensional Euclidean space (x) supplemented by two 
Grassmann dimensions (6, 6), which we collectively de- 
note by 6. Here, ft = z'ft — i : y {c, c) and the x-dependence 
only appears through the fields U, c, c and ft The first 
factor on the right-hand-side of d20l > is trivially an ele- 
ment of SU(N) and so is the superfield *V. 

Defining further the "super gauge transform" A v as 
in ©: A'Y = 'VA/V' 1 + ^-'Vd/V-\ it is straightfor- 
ward to show that the averaging weight in (flOb takes the 
remarkably simple form: 



5, ft] = I tn 

Jx,e 



S v/ [A u ]+S FP [A u ,c,c. 

with the notation 

£ = j d d xd9d8g 1/2 (9,9), 



(21) 



(22) 



Civita symbol. Our proof follows standard arguments 
01 80 and consists in identifying all local terms of mass 
dimension less or equal to four0 in the effective action 
r compatible with the symmetries of the theory. 

Let us first list the symmetries of the action ( 1241 
that are realized linearly. Apart from the trivial global 
SU(2) color symmetry and the isometries of the R 4 Eu- 
clidean space, one has a shift symmetry for the c field 
(c — > c + const.), the net ghost number conservation 
(c — > e ,e c, c — > e~' e c) and the isometries of the curved 
Grassmann space. The latter only impact the super- 
fields: r V k ^> r V k + X%<V k where X 6 - is one of the five 
independent Killing vectors on the Grassmann space 



3 This relies on Weinberg's theorem and assumes, in particular, that 
the free propagators decrease sufficiently fast at large momentum. We 
show in the next section that all free propagators behave as l/p 2 at 
large p, which is a sufficient condition. 
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lfl6Tl . At the level of the effective action, these sym- 
metries simply imply that terms involving Grassmann 
variables should be written in a covariant way: integrals 
always come with proper integration measure, see d22l i. 
and derivatives are contracted with proper tensors 11611 . 
An important remark to be made is that these transfor- 
mations apply to each individual replica superfield 'Vj, 
independently of the others. This implies that each such 
superfield always comes with its own set of Grassmann 
variables. There is also a discrete symmetry under per- 
mutation of the replicas: "V* <-> "V/, k, I — 2, . . . , n. 

These transformations are also symmetries of the ef- 
fective action T and directly constrain the possible di- 
vergent terms. We also use the fact that the linear term 
involving the field h is not renormalized: SF/Sh" - 
8S/6h a = id M A" and that the choice of the replica k = 1 
singled out in (lT8T>-(fT9]l being arbitrary, the divergences 
associated with the fields c, c and h are the same as those 
associated with q, q and hi for k > 2@ 

The action (1241 also admits non-linear symmetries. 
One is a BRST-like symmetry, corresponding to the in- 
finitesimal transformation: 



d,c a 



■ gO* AC 



sc a = -i°.f abc c b c c , 
2 

sV k = -ig ( {V k c, k = 2,. 



sc a = ih a , 

s(ih a ) = /3 Q c a , (25) 



It is useful to employ the decomposition of SU(2) ma- 
trices in terms of a unit 4-vector: 



"V* = n° k l + info" , 



(26) 



where (nf) + nfra? = 1 . The BRST transformation of 



the constrained superfield 'Vk then reads: 



,° = ^-nfc" and snl 



go 
2 



(-n° k c a + e abc n b k c c ) 



(27) 

Note that, just as "V*, n? and n a k are superfields, that is 
functions of (x, ff). In the following we choose n a k to be 
the unconstrained superfields and n? = (1 — n a k n a k y^ 2 . 

The action d24l is also invariant under 'V* — » U'Vk, 
where the SU(2) matrix H = U(ff) can be local in 
Grassmann variables. This symmetry is nonlinear since 
<Vi is a constrained superfield. Each replica super- 
field can be transformed independently of the others and 
there are thus 3 x (n - 1) such symmetries. The infinites- 
imal transformations read = iduffVi or, in terms 



of the decomposition 



8fo° = -dm" and S a k n h = 6 kl (<f fe n° + t^nf) . (28) 

To derive ST identities associated with the non- 
linear symmetries d25l and (1281 . one introduces (su- 
per)sources coupled to both the (super)fields and their 
variations under s and 6 a k . We define |f| 

J, = J^A a M +ri a c a +c a j 1 a +ih a M a +K^sA^+L a sc a } 

k=2 Jx 'i 

and consider the Legendre transform of the functional 
W — In J Dcf>e~ s+S ' with respect to the sources 7°, r/ a , 
ff, M a and P k , where <p collectively denotes the (su- 
perfields A"^, c a , c a , h" and n a k . It is a straightforward 
procedure to derive the desired identities. We obtain, 
for the BRST symmetry d25l-(l27l. 



I 



f 5Y 6T 6F 6T , „ 8Y . 6T 
< — — + ——-ih a — -foe" — 



[SA^SK" 5c a 5L a 5c a 



6ih a 



5Y 1 ST 6T 



+ S Xfl { 8 1,2 (i) ^ 5&1 + g(0) Snl 8QI ) ° 
and, for the symmetries 6% (no sum over k) 



(30) 



-n°— + O a — + e abc 



ST 



ST 



*5n' 



,1/2 



(6)1*4 + 



i sr si 



(31) 



where the factors g(ff) account for the curved Grass- 
mann directions. Observe that the identities (|3TT > are 
not integrated in 0, which follows from the fact that the 
symmetries 5 a k are local in Grassmann space. 

Dimensional analysis shows that the Grassmann vari- 
ables have mass dimension -1, which implies that J" . 
has dimension -2 in d = 4. This enables one to de- 
termine the dimension of the various fields and sources, 
summarized below together with their ghost numbers. 

It is then straightforward, although too lengthy to be 
reproduced here, to prove the renormalizability of our 
theory. One first writes the possibly divergent (local) 



4 For instance, upon the change of variables A — > A Ul , Uu — » 
UtU^ 1 for k > 2, JJi — > UZ 1 and c <-> C2, c <-» ci and h «-» hi, 
one gets that it is now the replica k = 2 which is singled out. 



Note that the symmetries * and S' k are not nilpotent. However, 
the transformations s 1 , s5 k = S a k s and 5^5* can be fully expressed in 
terms of the (super)fields at hand and their variations under either s or 
S" k . Therefore they do not require independent (super)sources. 
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field 


A 


c 


c 


ih 


n 


K 


L 


P 


G 


dim. 


1 


1 


1 


2 





2 


2 


2 


l 


ghost nb. 





1 


-1 








-1 


-2 





-l 



contribution to the effective action as Fdi v = f -C<Hv, 
where _£div contains all terms of mass dimension less or 
equal to four that are compatible with the linear symme- 
tries listed previously. Further imposing the constraints 
( f30b and OTI ). we get, after some calculations, 



div — -£>l ^u^r^u ^ ^ 



k=2 



Q>1 



Q>1 



(32) 



with the renormalized BRST variation: 



(33) 



and where the vanishing source terms can be written as 
£, = ^-{d,A a v - d v Al + goZ^Afaf + ^(A;f 



4Z,. 

d^fc'd^ + go^Afr) 



■ ih^Al 



and 



£ k = Z c tx\A 



goZ c 



(34) 



(35) 



We see that Fdiv has the very same structure as the 
original classical action up to the two (divergent) con- 
stants Za and Z c . It is remarkable that Tdiv only in- 
volves two divergent constants as in the usual FP Lan- 
dau gauge. Notice finally that each replica d35l > gives 
a gluon mass contribution fi()Z c (A a ^) 2 jT. after integration 
over Grassmann variables. Each such term is renormal- 
ized just as the one in (l3~4l . as expected from the replica 
symmetry. 



5. Feynman rules 

Perturbation theory is most conveniently formulated 
in the supersymmetric formalism, which makes trans- 
parent the (dramatic) consequences of the supersymme- 
tries - the isometries of the curved Grassmann space 
- for loop diagrams. To formulate Feynman rules, we 
parametrize the constrained SU(A0 superfields *Vj as 



r V k = expj/goA* 



(36) 



propagators of the theory, written below in Euclidean 
momentum space. Note that, because of the curvature 
of the Grassmann subspace, it is of no use to intro- 
duce Grassmann Fourier variables. The gluon propa- 
gator reads: 



[Al{p)A b v (-p)\ 



p 2 + nfio 



(37) 



where the square brackets represent averaging with the 
action (124-b . It is transverse in momentum, as a result of 
the Landau gauge condition, and massive, with (bare) 
square mass = n/?o, as a result of our particular 
gauge fixing procedure. Each replica contributes /?o to 
the square mass as already clear from Eq. ( fT9b . The 
ghost propagator assumes the usual form: 



[c a (p)-J(-p)] = 6 ab /p 2 . 
The superfield propagator is given by 

[A a k (p,0)Al(-p,0!)] = 6 ab 5 kl 5{0-i)lp 2 , 



(38) 



(39) 



where 6(6-0') = g- 1/2 (0)(0-O')(6-0')isthe covariant 
Dirac delta function on the curved Grassmann space: 

Jg8(i-£)f(@=f(2)B 

The vertices of the action d24l > which do not involve 
the superfields are the same as for the usual FP Lan- 
dau gauge. On top of the latter, there are vertices with 
an arbitrary number of legs and either zero or one 
gluon leg. Note that such vertices always involve the 
same replica. Coupling between different replicas can 
only arise through loop diagrams with gluons. 

The structure of the propagator ( 1391 leads to a dra- 
matic simplification of Feynman rules. Consider a loop 
of Ak superfields with p vertices insertions. It is easy 
to see that vertices involving A* superfields are local in 
Grassmann variables. Thus the loop involves the fol- 
lowing integral over Grassmann variables: 



6(6, -0,)-.. $(0-00 = 0. 



(40) 



We conclude that loops of At superfields vanish fl9h . 

This observation has two important consequences. 
First, the only source of dependence in the number n 
of replicas is the bare gluon mass in ( l37l l. Second, cor- 
relators or vertex functions involving only the fields A, 



where the superfields A^ are unconstrained. 

Expanding the action (l24l > in powers of the (su- 
perfields to quadratic order, we obtain the various free 



6 The correlators involving h are [h"(p)h b (-p)] 
[h"(p)A b II (-p)} = 6 ab Pll /p 2 and [h a (p)K b k (~p,0)} = -id' 
do not enter in loops since there is no vertex with a h leg. 



= 6"%/p 2 , 
ab /p 2 . These 
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c and c do not receive any contribution from replica su- 
perfields. They are thus obtained from the very same 
diagrams as in the FP Landau gauge, with usual YM 
vertices and with propagators given by ( l37T l and d38b . 
These are nothing but the Feynman rules of a massive 
extension of the FP gauge-fixed theory, which is a par- 
ticular case of the CF model. In particular, we recover 
the standard FP Landau gauge for /Jo = 0. 

It is interesting to consider the above result from a 
slightly different angle. Consider integrating out the su- 
perfields *Vk in (TT~8T > to obtain an effective theory for the 
YM and FP fields (A, c, c, h) which takes into account 
the average © over Gribov copies. Eq. d40b above im- 
plies that the corresponding functional integral receives 
no loop contribution and is thus (perturbatively) exact at 
tree levelQ Still the result of such an integration yields 
a highly non trivial effective theory: apart from a gluon 
mass contribution, it contains non local contributions of 
arbitrary order in the gluon field. However, it is easy to 
show that the corresponding vertices are actually longi- 
tudinal in momentum space with respect to at least two 
of their gluons legs. It follows that they do not con- 
tribute to standard YM correlators involving gluon or 
(anti)ghost fields since there they are to be connected to 
gluon propagators, which are transverse. Again, as long 
as one considers standard YM correlators, the only ef- 
fect of the non trivial gauge fixing procedure considered 
here is a mass term for the gluon, as in the CF model. 



6. Renormalization group flow 



The gluon and ghost two-point vertex functions have 



been computed at one-loop in the CF model 111 U 11211 
with dimensional regularization. The corresponding 
one-loop expressions for the present theory can thus be 
obtained from these papers by replacing the bare gluon 
mass mi — > n/3o. 

A non trivial issue concerns the order in which renor- 
malization and the limit n — * should be performed^ 
Due to the observation that the number n of replicas 
only appears through the bare gluon mass, the naive pre- 
scription to take first the n — > limit reduces to the 
usual FP gauge fixing. This is physically unsound since 
we expect ghost and gluon correlators to depend on /Jo. 



A possible renormalization scheme would be to intro- 
duce a renormalized parameter as fio = Zp /J. But again, 
this leads to an unsatisfactory n — > limit. Instead, 
we can make use of the fact that only the combination 
m? = nfio appears to implement the IR safe renormaliza- 
tion scheme proposed in fl2ll . We introduce the renor- 
malized fields and constants A = yfZ^Ar, c = ^Z c c r , 
c = ^JZZc r , go = Z g g and mi = Z m %m 2 and impose the 
following renormalization conditions 



rf(p=/i) = m z +// 



and 



Zs^Z~aZ c - 1, Z„{lZ^Z c - 1, 



(41) 



where T® and F?' denote the one-loop gluon and ghost 
renormalized two-point vertex functions. 

The renormalization group (RG) /J-functions for the 
coupling g and the mass m can be found in 1 12y. Here, 
we discuss the general properties of the RG flow. First, 
we find an ultraviolet (UV) attractive fixed point at m = 
and g = 0. This implies that upon removing the UV 
regulator (continuum limit) both the bare coupling go 
and the bare mass mo vanish. It is interesting to note that 
this behavior is compatible with taking the limit n — > 
at fixed Po along with renormalization^ Thus we may 
regard the various RG trajectories which emerge from 
this UV fixed point as corresponding to various choices 
of the gauge-fixing parameter /Jo. 

For one of these trajectories m = during the whole 
flow, which simply corresponds to the standard FP re- 
sult /Jo = 0. We show in Fig. Q] different trajectories ob- 
tained by integrating numerically the one-loop f3 func- 
tions. There are two qualitatively distinct families of 



60 



40 
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Landau / 
pole / / 
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\ 

m 

< \ 


J£^^~^ No Landau pole 

C— 3" 



2 4 ^6 8 10 

Figure 1: One-loop RG flow in the plane (m = m//i,g). 



7 That this is the case can be understood by remembering the fact 
that this integral is a way to sum over the solutions (/, of the classical 
field equation d^Ajj = at fixed A, see Eq. UOt . 

8 A similar issue arises in the context of disordered systems with 
the thermodynamic limit. The common understanding is that the latter 
should be performed first, see e.g. Ilal . 



trajectories divided a separatrix (dashed line on Fig. fT}. 
The trajectories above the latter, which include the stan- 
dard FP case, are characterized by the presence of a 



9 A similar observation was made in a different context in [20]. 
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Landau pole as one run towards the IR: the coupling g 
diverges at a finite momentum scale /y. In constrast, for 
trajectories below the separatrix, the coupling constant 
remains bounded all the way to fi — > 0. These are par- 
ticularly interesting because the effective coupling con- 
stants are well-defined at all scales which means that 
the IR properties of the theory may be accessible by 



perturbation theory. It was actually shown in 111 ll 11211 



that one of these trajectories leads to two-point ghost 
and gluon correlation functions in rather good agree- 
ment with lattice simulations. This might seem surpris- 
ing at first sight since the gauge fixings are different. 
However, for /3q not too small the saddle points are sup- 
pressed in Eq. ([2} and for/Jo not too large all minima are 
equiprobable. Therefore, if these conditions have some 
overlap for some range of f3o, we expect the usual lattice 
implementation of the Landau gauge to be similar to the 
gauge fixing proposed here. 

The separatrix ends in an IR fixed point where the 
coupling constant is large. Its location should thus not 
be taken seriously. However, we expect the qualitative 
features of the flow described above to be stable against 
higher order corrections in perturbation theory. Indeed, 
on the one hand it is expected that the existence of a 
Landau pole for the trajectory m — is a property valid 
at all orders and, on the other hand, trajectories lying 
close to the g — axis should not be drastically influ- 
enced by higher order corrections. 

7. Conclusion 

We have devised a well-motivated, first-principle for- 
mulation of YM theories in Landau gauge, which takes 
into account the issue of Gribov copies in a consistent 
way and avoids the Neuberger zero problem. Perturba- 
tion theory shows no Landau pole and is thus potentially 
under control at all momentum scales. This questions 
the standard understanding that the low-momentum sec- 
tor of YM theories is genuinely nonperturbative and 
opens the exciting possibility that the highly nontriv- 
ial IR physics may be accessible by perturbative meth- 
ods. This opens the way to further investigations includ- 
ing e.g. higher-order corrections, finite temperature, or 
quark dynamics. 

Another line of development concerns the study of 
classical solutions of the gauged-fixed theory proposed 
here. It was pointed out long ago that the classical equa- 
tions of motion of a similar gauged nonlinear sigma 
model present non-trivial vortex solutions which may 
be related to confinement [21]. 

Finally, it is an interesting question whether the gauge 
fixing procedure proposed here, based on averaging 



over Gribov copies rather than choosing a single one, 
can be (approximately) implemented in numerical sim- 
ulations on the lattice, where a study of the influence of 
the saddle points would be necessary, or in nonpertur- 
bative continuum approaches such as Schwinger-Dyson 
equations or the functional RG [22]. 

We thank G. Tarjus and N. Wschebor for fruitful in- 
teractions during the early stages of this work and for 
valuable suggestions concerning the limit n — > as well 
as L. Cugliandolo, Vi. Dotsenko and M. Mezard for 
useful discussions. 
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